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Abstract 

In the present paper we describe the dynamics of the revised rigid 
body, the dynamics of the rigid body with distributed delays and the 
dynamics of the fractional rigid body. We analyze the stationary states 

■ for given values of the rigid body's parameters. 
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. 1 Introduction 

■ In mechanical problems, the dynamics of the rigid body has an important 
role. In many papers M. Puta analyzed the dynamics of the rigid body with 
control and he obtained important results. 

Recently, the dynamics of revised rigid body, the dynamics of the rigid 
body with distributed delay and the dynamics of the rigid body with frac- 
tional derivative have been studied. The last two aspects represent the dy- 
namics of the rigid body with memory. 

Our paper studies the dynamics of the revised rigid body obtained by 
a metriplectic structure which is canonically associated. We define the dy- 
namics of the rigid body with distributed delays. For the Euler-Poincare 
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dynamics of the rigid body we analyze the hnearized system in an equihb- 
rium point. We obtain the the existence conditions for the Hopf bifurcation 
with respect to the parameter of the repartition density which defines the 
distributed delay. Also, we define the fractional rigid body dynamics using 
the Caputo fractional derivative. 

2 The revised differential equations for the 
rigid body 

The differential equations for the rigid body in W are described by a 2- 
antisymmetric tensor field P and the Hamiltonian function h given by: 

/ \ 

P{x)^{P'^{x))^{ -x^ x' ], h^-{a^{x'y+a2{x'f+as{x'f), 
\ x^ -x' J ^ 

(1) 

where {x^,x'^, x^Y ^ ^ ^ = 1) 2, 3, Oi > 02 > a^. These differential 

equations are: 

x{t)^P{x{t))V,h{x{t)), (2) 

where x{t) = {x^{t),x^{t),x^{t)Y and Vxh is the gradient of h with respect to 
the canonical metric on R^. The differential equations (2) have been studied 
by M. Puta in [8]. 

Let M? be the space interpreted as the space of body angular velocities 
VL equipped with the cross product as the Lie bracket. On this space, we 

consider the standard Lagrangian kinetic energy L{VL) = -I ■ Q, where I = 

diag{Ii, I2, 13) is the moment of inertial tensor, so that the general Euler- 
Poincare equations become the standard rigid body equations for a freely 
spinning rigid body: 

m^{i-n)x n. (3) 

If M = / • Q is the angular momentum, then, (3) is: 

M = M X O. (4) 
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If M = {Iix{t),l2yit),l3z{t)f, Q = {x{t),y{t),z{t)f from (4) results: 
x{t) = ^-^^y{t)z{t), y{t) = k^^^t)z{t), z{t) = ll^x{t)z{t), 

-tl -'2 

(5) 

with Ii > I2 > Is- The differential equations (5) have been studied by M. 
Puta in [9]. The revised differential equations for the rigid body given by 
(2) have been studied in [5]. They are described by P, h and the tensor 
fields g — {g'^^{x)), where P and h are given by (1) and g is defined by 
g{x) = {g'^x)), 

and the Casimir function of structure Poisson P is: 

c{x)^\{{x'f + {x'f + {x'f). 

The structure , g,h,c) is called a metriplectic manifold of second 
kind. The revised differential system associated to (2) is given by: 

x{t) - P{x)V:,h{x) + g{x)V:,c{x) . (6) 

Proposition 2.1. (i) The differential equations (2) are given by: 



x\t) = {a2-a3)x^{t)x^{t),x^{t) = {a3-ai)x\t)x^{t),x^{t) = {ai-a2)x\t)x'^{t); 

(7) 

(ii) The differential equations (6) are given by: 

x^{t) = {a2-a3)x'^{t)x^{t)+a2{ai-a2)x^{t){x'^{t))'^ + a3{ai-a^)x^{t){x^'{t)Y 
x'^{t) = {a3-ai)x^{t)x^{t) + a3{a2-a3)x'^{t){x^{t)f + ai{a2-ai)x'^{t){x\t)f 
X^{t) = {ai-a2)x^{t)x'^{t) + ai{a3-ai)x^{t){x^{t)f + a2{a3-a2)x^{t){x'^{t)f; 

(8) 

The equilibrium points of the system (7) are studied in [8], and the equi- 
librium points of the system (8), in [5]. 
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3 The differential system with distributed de- 
lay for the rigid body 

We consider the space M^, the product x = {{x, x),x G M^, x G M^} 
and the canonical projections TTj : x ^ M"^, z = 1, 2. A vector field 
X e X(R^ X M^) satisfying the condition X{nlf) = 0, for any / G C7°°(M3) 
is given by: 

d 



X{x,x) = Xi{x,x) — . 
The differential system associated to X is given by: 



x'{t) = X\S;{t),x{t)),i = 1,2,3. 

A differential system with distributed delay is a differential system associates 
to a vector field X G X(R''^ x R^) for which X{Tr*J) = 0, for any / G C= 
and it is given by (1), where x{t) is: 

x{t) = / k{s)x{t — s)ds, 







where k{s) is a density of repartition. In what follows, we will consider the 
case of the following densities of repartition: 

(i) the uniform density with: 

{0, < s < a 
-, a < s < a + T , 
0, s > a + T 

where a > 0, r > are given numbers; 

(ii) the exponential density, with ka{s) = ae~"^, a > 0; 

(iii) the Erlang density, with ka{s) = a^se""*, a > 0; 

(iv) the Dirac density, with kr{s) = 6{s — r), r > 0. 
The initial condition is: x{s) = (p{s), s G (— oo,0], 

with (p : (— oo, 0] — > a smooth map. Some systems of differential equations 
with distributed delay in R^ were studied in [1], [2]. For such a system, we 
consider relevant the geometric properties of the vector field which defines the 
system, for example first integrals (constant of the motion), Morse functions. 
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almost metriplectic structure, etc. The differential equations with distributed 
delay for rigid body are generated by a 2-antisymmctric tensor field P on x 
that satisfies the following relations: P{nlfi, 112/2) = 0, P{Txlfi, 1^1/2) = 0, 
for all /i, /a G ^'^(M^) and h e C°^{R^ x M^). xhe differential equation is 
given by: 

X — P{x,x)V xh{x,x). (9) 
Let P(x, x) be the tensor field with the components given by: 






x^ 


—x^ 


—x^ 





x^ 


x" 


-x^ 






{P'\x,x)) = I x^ \ (10) 

and 

h [x, x) = ttix^x^ + a2X^x^ + a^x^x^ . (11) 
The differential equations (9) are: 

x^{t) = a2X^{t)x^{t) - a^x^{t)x^{t), 
x^{t) = a^x^{t)x^{t) - aix^{t)x^{t), 
x^{t) = aix^{t)x^{t) - a2X^{t)x^{t). 

The differential system: 

X — P{x, x)'Vxh{x, x) + g{x, x)'Vxc{x, x), (12) 

where the components of g{x, x) are: 

dh{x,x)dc{x,x) ^ dh{x x)dc{x x) 

(13) 

is called the revised differential system with distributed delay associated to 
the differential system (9). 

Prom (10), (11), (13) and c{x,x) = l{x^f + x'^x'^ + l{x^f results: 

/ — a|a;^x^ — a^x^x^ aiOaX^a;^ aia^x^x^ \ 

{9^'')— aiCax^T^ —a\x^x^ — a2,x^x^ a2a^x^x^ 

^ a\a^x^x^ a2azxP'x^ —a\x^x^ — a\x'^x^ J 
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The differential system (12) is given by: 

x\t) = a2X^it)x^t) - a3X^{t)x^{t) + aia2x\t){x^{t)y, 

x''{t) = a3x\t)x^{t) - aix\t)x^{t) - alx\t)x\t)x''(t) - alx''(t)x^t)x^t), 
x^lt) = aix\t)x'^lt) - a2x\t)x\t) + a2a3X^{t)x^x^ . 

(14) 

For the Dirac distribution, system (14) was analyzed in [6]. The other 
types of densities will be analyzed in our future papers. 

The Euler-Poincare equation for the free rigid body with distributed delay 
is defined by: 

M = M xn + aM X {M xn) 

where M = (/ix(t),/2i/(t),/32;(t))^,l] = (x(t),y(t),^(t))^,0= (x(t),y(t),5(t))^, 
M = m, h > 0, /2 > 0, /g > and a G M. 

The equilibrium points of our system are Qi = {j^, 0, 0)-^, fl2 = (0, ^, 0)-^, 
^3 = (0,0,|^)^, m e R*. 

Proposition 3.1. The equilibrium point fli has the following behavior: 
(i) The corresponding linear system is given by: 

U{t) = AU{t) + aBU{t) 

where U{t) = {u^{t),u^{t),u^{t))^ and 

/ 0\ /OO 

A=i ],B= [ 

\ ^hTT^ / V ^7777^^ 
(a) The characteristic equation is: 

-'2 J^3 -'l-'2-t3 

(Jl-/2)(/3-/l) 



-m'] = 0; 



/f/2/3 

(Hi) On the tangent space at Qi to the sphere of radius the linear 
operator given by the linearized vector field has the characteristic equation: 

12 3 2 3 

(/l-/2)(/3-/l) 
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(iv) If Ii > I2, h > h and k^'^^X) = e~'^^, t > 0, for < t < Tc, where 

^ h{h{h-h)+l2{h-h)) 

?>\a\m?{h-h){h-h) ' 
then the equilibrium point Qi is asymptotically stable. 

The analysis of the equihbrium point fti for the Dirac density is given in 

[1]- 

4 Fractional differential systems for the rigid 
body. 

Generally speaking, the fractional derivative, Riemann-Liouvillc frac- 
tional derivative and Caputo's fractional derivative are mostly used. In the 
present paper we discuss the Caputo derivative: 

D^x{t) ^ I^-'^ij^rxit), a>0, 

where m — 1 < a < m, (-7-)'" = -7- ° •■• ° -7-, -^^ is the /?th order Riemann- 

at at at 
Liouville integral operator, which is expressed as follows: 

I^x{t)^:^j\t-sf-'x{s)ds, /3>0. 

In this paper, we suppose that a e (0, 1). 

Examples of the fractional differential systems arc: the fractional order 
of Chua's system, the fractional order of Rossler's system and the fractional 
Duffing oscillator. The geometrical and mechanical interpretation of the 
fractional derivative is given in [7]. The geometry of fractional osculator 
bundle of higher order was made in [3] using the fractional differential forms 
[4]. 

A fractional system of differential equations with distributed delay in 
is given by: 

Dtx{t) = X{x{t),x{t)), ae(0,l), (15) 

where x{t) = {x^{t),x'^{t),x^{t))'^ e R^. The linearized of (15) in the equi- 
librium point xo , {X{xq,xo) — 0) is given by the following linear fractional 
differential system: 

D^u{t) = Au{t) + Bu{t), (16) 
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where A = (U S = (If 

The characteristic equation of (16) is: 

A(A) = detiX^I -A- k^'\\)B), (17) 

where k^^\X) = k{s)e^'^^ds. 
Prom (16) we have: 

Proposition 4.1. ([4]) (i) If all the roots of characteristic equation A(A) = 
have negative real parts, then the equilibrium point xq of (15) is asymptotically 
stable; 

(a) If k{s) is the Dirac distribution, the characteristic equation (17) is 
given by: 

A(A) = rfet(A"/ -A- e-^^B) = 0. 

If T = 0, a E (0, 1) and all the roots of the equation det{\I — A — B) = 
satisfy \arg{X)\ > then the equilibrium point xq is asymptotically stable; 

(Hi) If a e (0.5, 1) and the equation det{XI—A—Be~^'^) — has no purely 
imaginary roots for any t > 0, then the equilibrium point Xq is asymptotically 
stable. 

For the following delayed fractional equation (see [4]) 

D^x{t) = ax{t - r) (18) 

where a G (0, 1), a G M and r > the stability condition is: 

lfa<0, (-a)^ 7^ i((2A; + l)7r - fTr) and (-a)^ -^((2A; + l)7r - f tt), 

k E'Z, then the zero solution of (18) is asymptotically stable. 
For the following delayed fractional equation (see [4]) 



D^x{t) = y{t) - k,x{t) 

Dty{t) = -{kr + k2)y{t) + x{t-T), 



(19) 



where a G (0, 1), /ci > 0, /c2 > 0, r > 0, the stability condition is: 

If ki > 0, k2 > — k, then the zero solution of system (19) is asymptot- 
ically stable. 

For / G C~(M3), by D° /, D° /, D'^^f we denote the Caputo partial 
derivatives defined by: 

(20) 
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where (x*) are the coordinate functions on IR, and (^), i — 1,2,3 is the 



^ dx 

canonical base of the vector field on R". 
Prom (20) results: 

^A^) r(l + 7-a) '^-'^"'^ 

Let X"(M'^) be the module of the fractional vector fields generated by 
the operators {-D"i,i = 1,2,3} and the module !D(]R^) generated by 1-forms 
{rf(x^)", i = 1, 2, 3}. The fractional exterior derivative rf" : C°°(R=^) ©(R^) 
is defined by: 

d'^{f) = d{xrD^.{f). 

Let P G X"(R"'') X X"(R^) be a fractional 2-skew-symmetric tensor field 
and rf"/, d'^g G D(R3). The bilinear map [•,•]" : C°°{R^) x C°°(R3) 
C°°(R3) defined by: 

oo /'TroSN 



[f,gr = Pid''f,d''g),yf,g G C^{R') (21) 

is called the fractional Leibnitz bracket. 

If P = p'W^i ® D^j then, from (21) it follows that: 

[f,gr = P'W:.fD:,g. 
From the properties of the fractional Caputo, results: 

ifh^ar = E ( ^ ) p''iD7'f)iD:^9) {J-jh 

a 

If P is skew-symmetric we say that (R'^, [•,•]") is a fractional almost Pois- 
son manifold. If a 1 then we obtain the concepts from [8]. 

For h G C°°(R^), the fractional almost Poisson dynamic system is given 

by: 

D-^x\t) = [x'{t),h{t)]'', where[x\ /i]" = P'W^.h. (22) 
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Let P be a 2-skew-symmetric fractional tensor field and 9 a 2-symmetric 
fractional tensor field on R^. We define the bracket [•,•]" : C°°(R^) x 

C°°(M3) ^ (:7°°(M3) by: 

a a 

The structure (M, P,9, [■,■]) is called fractional almost metriplectic man- 
ifold. The fractional dynamic system associated to /i e C°°(M^) is 

D^x'it) = [x'{t), hit)]'', where[x\ h]'' = p'W^^h + g'W^^h. (23) 
If we define the bracket [•,(•,•)]" : C°°(M3) xC~(M3) xC°°(M3) ^ C^^R^) 

by: 

[/, ih, hT = P{d"f, d'^h) + gid'^f, d''h2)yf, h, h e C^iR'), 

a 

then, the fractional vector field Xh^hz defined by: 

is called the fractional almost Leibnitz vector field associated to the functions 
The fractional almost Leibnitz dynamical system is given by: 

D^x'it) = pW^,hi + g'W^M- 

a . . a 

Let P = (P*-'), 9 ~ (gij) be the fractional 2-tensor fields on K'* and 
h e C°°(M3) given by: 

hi = f(^^K(^')"^' + <^2{x'r^' + asix'r^']. 

Proposition 4.2. (i) The fractional dynamic system (22) is: 

D^x^ = (as - as)x'^x^, Dfx"^ = (ag - al)a;^x^ D'^x^ = (ai - a2)x^x'^; (24) 

(a) The fractional dynamic system (23) is: 

Dfx^ = (a2 — a3)x'^x^ + a2(ai — a2)x^{x'^Y + 03(01 — a3)x^{x^Y 

D^x^ = (03 — ai)x^x^ + 03(02 — 03)x^(a:^)^ + 01(02 — Oi)a;^(x"^)^ (25) 

D^x^ = (oi - a2)x^x'^ + 01(03 - ai)a:^(a;^)^ + 02(03 - 02)x^(a:^)^; 

10 



(in) The fractional dynamic systems (24) and (25) have the equilibrium 
points Mi(m, 0, 0), M2(0, m, 0), M^{Q, 0, m), m E R* ; 
(iv) The characteristic equations for (24) are: 
in Mi(m,0,0); A°(A2° + (ai - a3)(ai - a2)m'^) = 0, 
m M2(0,m,0).- A"(A2° - (ai - asKaa - ) = 0, 

in M^lO , , m) : A"(A^" + (ai — a3)(a2 — a3)m^) =0; (^w^ The characteristic 
equations in Mi (m, 0,0), M2(0,m,0), M3(0,0,m) for (25) are: 

A"(A^° — ai(a2 + 03 — 2ai)m^A" + (ai — a3)(ai — a2)m^{a\m^ + 1)) = 0, 

A"(A^° — 02(01 + 03 — 2a2)m^A° — (oi — a2)(a2 — a'i)m^{a\m^ + 1)) = 0, 

A"(A^° — 03(01 + 02 — 2a3)m^A" + (ai — 03) (02 — a^)rin?' {a^rn^ + 1)) = 0; 

The above findings allow the analysis of the equilibrium points with re- 
spect to the parameters of the characteristic equations. 

Equations (24) are called the fractional equations of the rigid body and 
equations (25) are called revised fractional equations of the rigid body. If 
o; — > 1, the results from Proposition 4.2 lead to results from [8]. 

For tti := 3, 02 := 2, 03 := 1, a = 1 the dynamics [x^ {t) , {t) , {t)) of 
(24) is given in figure Fig.l and for a = 0.82 in figure Fig. 2. The numerical 
algorithm used is Adam-Moulton-Bashford. 

Fig.1.alfa=1 Fig.2.alfa=0.82 




5 Conclusions 

In the present paper we present the dynamics of the rigid body with 
memory. The memory was described by the variables with distributed delay 
and by the Caputo fractional derivative. 
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